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TECHNICAL MEMORANDIM 1^79 

A METHOD OF QUADRATURE FOR CALCULATION OF THE LAMINAR 
AND TURBULENT BOUNDARY LAYER IN CASE OF PLANE 
AND ROTATIONALLY SYMMETRICAL FLOW* 

By E. Truckenbrodt 


1. SlMdARY 


For calculation of the characteristic paraineters of the boundary 
layer (momentum- loss thickness and form parameter for the velocity pro- 
file), two quadratiire formulas are given which are valid for the lam i nar 
as well as for the turbulent state of flow. These formulas cover both 
the two-dimensional and the rotationally symmetrical case . 

The calculation of the moraentimi-loss thickness is carried out by a 
simple integration of the energy theorem. The equation for the form 
parameter is obtained by coupling of the momentum theorem with the energy 
theorem. Knowledge of the derivatives of the velocity distribution and 
of the radius of the body along the length x is not necessary. 


2. INTRODUCTION 


The calculation of the laminar and turbulent boundary layer in case 
of pressure drop and pressure rise is the decisive problem in the deter- 
mination of the flow loss in channels and pipes and of the flow drag of 
bodies . We shall treat below the plane problem as well as the rotationally 
symmetrical one . In these considerations , we shall limit ourselves to 
incompressible flows . 

Whereas theoretical treatment of the laminar boundary layer has 
been fundamentally clarified, one is still dependent on semiempirical 
connections for the treatment of the turbulent boundary layer. Since the 
exact methods in the laminar case are rather troublesome, various approxi- 
mation methods have been developed just like for the turbulent case . All 
these methods are based on the momentum equation given for the first time 
by Th. von K^m^ (ref. 1) . 

*"Ein Quadraturverfahren zur Berechnving der laminaren und turbulenten 
Reibungsschicht bei ebener und rotationssymraetrischer Stromiing." Ingenieiu*- 
Archiv, Band XX, Viertes Heft, 1952, pp. l6-228. 
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The approximation method for the plane laminar case has been per- 
fected by K. Pohlhausen (ref. 2). Later on, it was essentially improved 
by H. Holstein and T. Bohlen (ref. 3) : they introduced as the desired 

parameter not the boundary- layer thickness as Pohlhausen had done but 
the momentimi-loss thickness. More recently, K. Wieghardt (ref. 4) der- 
ived an energy equation in addition to the momentvim equation which he 
uses for developing a two-parameter calculation method . A . Walz (ref . 5) 
simplified this latter method by reverting to the one -parameter condition 
as in the methods of Pohlhausen and Holstein -Bohlen. The rotationally 
symmetrical method which is analogous to the Pohlhausen method was indi- 
cated by S. Tomotika (ref. 6) . Its simplification in the sense of the 
Holstein-Bohlen statement was carried out by F. W. Scholkemeyer (ref. J) . 

A method for calculation of the plane turbulent boundary layer was 
indicated for the first time by E. Gruschwitz (ref. 8) . Aside from the 
momentum theorem, Gruschwitz also uses a semiempirical equation obtained 
from certain energy considerations for determination of a form parameter 
marking the velocity profile in the boundary layer; this form parameter 
characterizes the sensitivity to separation of the boundary layer. The 
empirical parameters appearing in Gruschwitz' were investigated once more 
by A. Kehl (ref. 9) • Another method which was similar to Gruschwitz' 
method was developed by E. Buri (ref. 10). In the United States, a 
method by A. E. von Doenhoff and N. Tetervin (ref. 11) proved to be use- 
ful. In this method, too, an empirical equation, in addition to the 
momentum theorem, is used for determination of a form parameter char- 
acterizing the velocity profile. Starting from this report, H. C. Garner 
(ref. 12), England, developed a method which is superior to that of 
V. Doenhoff -Tetervin in its numerical evaluation. The transfer of the 
momentum theorem to the rotationally symmetrical case was performed by 
C. B. Millikan (ref. 15). 

Owing to recent investigations by H. Ludwieg and W. Tillmann (ref. l4) 
and J. Rotta (ref. I 5 ) oh the theoretical properties of turbulent flows, 
particularly of the wall shear stress and the energy loss in the boundary 
layer, it is possible to find a better basis for and to improve the 
existing semiempirical methods . 

It is the aim of the present report to develop a calculation method 
which is equally valid for the four cases of laminar and turbulent as 
well as plane and rotationally symmetrical flow. The most recent results 
will be taken into consideration. 


3 . THE MOMENTUM THEOREM AND ENERGY THEOREM 

As is well known, Prandtl's bo vindary- layer equations and Bernovilli's 
equation represent the fundamental equations for calculation of boundary 
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layers • As stated in the Introduction, we are going to consider incom- 
pressible flows. Since, however, the following derivations can be given 
for variable density p without particular difficulties, we shall make 
the transfer to incbmpressible flow only in the final result. 

We assime that, in the rotationally symmetrical case, the radius 
of the body or pipe R is large compared to the thickness of the boiindary 
layer &. The equations mentioned then read 



^yy by 


( 1 ) 


^(puR) + ^(pvR) ^ 0 

bx by 


(2) 


^x 


PaU 


dx 


(3) 


Therein, u and v signify the velocity components within the boundary 
layer 


0< y*i 5 


in X and y direction, U the velocity outside of the boundary layer 

(y ^ B) > figure 1. p is the pressxare, assumed constant across the 
bovindary-layer thickness, t the shear stress, and p the density of 
the flowing mediiam. For the plane (two-dimensional) case, the radius R 
is to be omitted in equation (2) as well as in all following formulas. 


The pertaining boundary conditions are 


II 

0 

11 

0 

II 

> 

T = Tq, 

P = PQ_ 

y = 6 

: u = U, 


T = 0, 

P = Pa 


If Vq 0, one is dealing with the case of suction (vq < 0) or of 
blowing (vq > 0) . 
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We now combine the three equations (l) to (5) adding to equa- 
tion (1) - -which has been m\iltiplied by - the eqxaation (2) which has 

been mviltiplied by u^‘*’^/(r + l)R; one may choose arbitrarily 
r = 0, 1, 2. . . 


One then obtains 


- ^(pu^+2r) + A(pur+lv) = (r + l)ur/paU (5) 

R dx oy y dx dyj 


The continuity equation (2) is integrated over the distance y from the 
wall, and one obtains 


pv 


= . i A 

R dx Jq 


puR dy + pqVq 


( 6 ) 


If we now assimie the velocity distribution u(y) and the shear -stress 
distribution r(y) in y direction to be known, we may integrate equa- 
tion (5) over y from y = 0 to y = 6. We then obtain, if we, more- 
over, substitute equation (6) 


1 d Tir+2Rf W & PO 2!0 _ _ 


p„iF‘^% 


(7) 


Therein, the newly introduced abbre-vLations signify 


dy 


T6 ^ 

r . 

.r-fli 

P ^ 

1 ^ 1 

~ ^ Pa U 

\uj 

= - (r + 1) 

n6 

/ S 

JL - 

Jq u 

Pa 


r-1 


dy > 


6r = 


- (r + 1) 



ul"" d, 


U/ ^y\PaU^ 


dy 


( 8 ) 
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Equation (7) is valid quite generally for laminar and tvarbulent, incom- 
pressible and compressible, plane and rotationally symmetrical flow with 
and without suction or blowing. As said before, the radius R is to be 
omitted in the plane case . 

If r is put equal to zero, there results the well-known momentum 
equation of von KirmAn. K. Wieghardt (ref. k) and J. Rotta (ref. 15) have 
shown that, for r = 1, equation (7) niay be given a physical interpreta- 
tion, namely that of the energy equation. We shall now assume that the 
flowing medium is incompressible, p = Constant, and also that neither 
suction nor blowing occur vq = 0 . Then r becomes 

r = 0: momentum theorem (v. KArmin, Pohlhausen) 


J_ ^u2r(|) , 61 = SL, 

U^R ^ U dx py2 


(9) 


with 


d = fo = 


= go = 





as momentum- loss thickness (lO) 


as displacement thickness (ll) 


= en as wall shear stress 

pu® ” 


r = 1: energy theorem (Wieghardt, Rotta) 


^ i^(u5R6) = 2 
^ pu5 


( 12 ) 
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with 


8 = 


fl = 



as energy-loss thickness (13) 


d + t 

pU^ 




T ^ Tu' 

pl^ SyVU 


4 y 


as shear-stress work 


1 


( 14 ) 


In the laminar case , the shear-stress work equials the energy converted 
into heat (dissipation d) . In case of tijrbulent flows, not all energy 
is converted into heat; one part still remains as .turbulence energy (t) 
which is, however, usmlly negligibly small. The fact that, in the case 
r = 1 , the second term on the left side of equation (7) vanishes is 
important since, with the density assumed to be constant, gj^ = 0 . 

If we now introduce the boundary -layer thickness ratios 



and H = S. 


we may write for equations (9) and (l 2 ) 


(15) 


— — (U^M) + H - 
u 2 r dx U 




(16) 


JL(u 5 i®^) = 2 ^ (17) 

U^R pU5 


We subtract equations (16) from (17) and find^ 

^his equation is obtained by partial integration from ( 8 ) . 

p 

For the laminar case, a corresponding formula has already been 
given by A. Walz (ref. 5 ) • la the turbulent case, too, one can show that 
the empirical equations indicated by E. Gruschwitz, A. E. von Doenhoff 
and N. Tetervin, and also H. C. Garner for calc \ilat ion of the form param- 
eter of the velocity profile (in Gruschwitz' paper q = 1 - (ud/u)^, in 
those of the others H) have a structure very similar to that of 
equation (18) . 
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dx U dx 


d. + t 


.U5 



( 18 ) 


As will be shown in the following section, the shear stress Tq^pU^ and 

the shear-stress work (d + t) /(i^ may be expressed as functions of the 
dimensionless quantities U'fl/v (Reynold.s number fomed with the momen- 
tijm thickness) and H. Furthermore, a one-parameter condition is valid 
in good approximation for the velocity distributions in the boundary 
layer; that is, there exists a fixed relation between the boundary-layer 
thickness ratios fl and H. With consideration of these facts, the 
equations (17) and (l8) represent two equations for the two ^lnknowns 
and H. 


In contrast to the existing methods, we shall calculate the momentum- 
loss thickness not from the momentum theorem (l6) but from the energy 
theorem (17) • This offers decisive advantages for the performance of the 
integrations indicated in section 4 . We resort to the momentum theorem 
in connection with the energy theorem for calculating the parameter H, 
equation (I8), which is characteristic for the velocity profile. 


4 . SHEAR STRESS, SHEAR-STRESS WORK, AND BOUNDARY-LAIER 

THICKNESS RATIO 


In order to be able to operate with the equations (17) and (I8), one 
must know the dependence of the shear stress and the shear-stress work on 
the quantities Ud/v and H, and likewise the connection H(H) . 


(a) Laminar Flow 

Assuming single -parameter velocity profiles we may write 



idience then follows 


Tq _ a2(H) 

pu2 



with 


(19) 
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and 


pu5" Hi 



( 20 ) 


H = 





( 21 ) 


For evalimtion of these formxalas, we select the so-called Hartree pro- 
files; these are the exact profiles for the velocity variation U ~ xm. 
In figure 2 , the quantities a, 3, H are plotted against H It is 
of interest that 3 is almost independent of H. For the case without 
pressure gradient, one has a = 0.220, 3 = 0.175> H = 2.6o. 


(h) Turbiilent Flow 

Recently, H. Ludwieg and W. Tillmann (ref. l 4 ) , as well as J. Rotta 
(ref. 15), dealt in detail with the determination of the wall-shear stress 
in case of turbulent boundary layers with pressure gradient . 

Ludwieg-Tillmann indicate for the range of the Reynolds numbers 
1 X 105 < U^/v < 4 X 10 ^ the following interpolation formula 


J[0 = 10-0 .678H (22) 

Ivi 


Rotta finds for the wall-shear stress the relation 



took over the numerical values of A. Walz (ref. 5 )* Ihe fol- 
lowing identities are valid: = e*, 3^ = D* and S= H^2* 
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*.c 


Therein 



signifies the shear velocity. 


Thus 



(24) 


is valid. In equation (25) 1 /k = 2.5 and B = B(Ii) is a function 

of the quantity It = — . The function B was evaluated and 

H u* 

graphically represented by Rotta. 

We calculated the shear-stress values for various values of H and 
plotted them against the Reynolds number Ud/v in figure 5. The agree- 
ment between the values according to Ludwieg-Tillmann and to Rotta is 
quite satisfactory. 


J . Rotta also dealt with the calctilation of the shear -stress work 
(dissipation and turbxilence energy) . For the dissipation, there applies 



f-)^ 

i Infn + G 

U/ 



(25) 


The fxuiction appearing therein, G = G(Ii), has been evaluated by Rotta. 
With consideration of equation (25), one may write for equation (25) also 



(25a) 


Since = f and H 1 JL^ the dissipation may also be 

represented as a function of the Reynolds number U^/v and of the form 
parameter H; this has been done in figvire 4. It is found that the 
differences for different values H are only slight; whereas, for the 
wall-shear stress, they are very considerable. We see from figvire 4 
that the dissipation may be approximated by the statement 

pu5 


(26) 
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A suitable value for n is n = l/ 6 . The dependence of the value 3 ^ 
on H is only slight so that we may assume for it the constant value 

3 = 0.56 X 10 “ 2 . 

For calculation of the turbulence energy, J • Rotta indicates an 
approximation formxila. With the aid of this formula, one can show that 
the turbulence energy is negligibly small compared to the dissipation; 
this fact has already been pointed out by Rotta- We put, therefore, 
for our f\irther calculations 

-^ = 0 (27) 

pu5 

The bomdary-layer thickness ratio H (energy-loss 
thickness/raomentum-loss thickness) may also be determined from J. Rotta 's 
results . The equation 


H = 


5 - 



(H - l)g 
H 


( 28 ) 


is valid. Therein, I 2 is a form parameter which is in a fixed relation 
(indicated oy Rotta) with We calculated according to equation (28) 

the values 3 in dependence on the quantity H and on the Reynolds num- 
ber W/v. The influence of the Reynolds number was shown to be vanish- 
ingly small. The result of our calculation is represented in figure 5* 

K. Wieghardt (ref. 16 ) also has dealt with the connection H(H) . He finds 


ii = — with A = 1.269 and B = 0-579 ( 29) 

H - B 


The resulting curve also has been plotted in figure Except for 
the larger values of H, the agreement with Rotta 's curves is satisfactory. 


5 • THE APPROXIMATION METHOD 


We shall simimarize the result of the previous section: The repre- 

sentations 



pU^ 



lit = ^ 
pu5 pu3 VI 


and H = h(H) 



NACA TM 1579 


11 


given in figures 2 to 5^ are valid. In the laminar case, the 
values Tq/pI^ and d/pU^ are inversely proportional to the Reynolds 
number • 


(a) The Calculation of the Momentum-Loss Thickness 

As mentioned above, we shall determine the momentum-loss thickness 
from the energy theorem (l?) • 

We write for the shear-stress work according to equations (20), (26), 

smd ( 27 ) 


<i -t- t _ 3 (h) 

pu5 / 


( 30 ) 


with n = 1 being valid for the laminar, n = 1/6 for the turbulent 
case . 


We now substitute this expression into equation (ij) and obtain with 


n 




(51) 


-^ = 2(1 + n)3(H)S”u5+2nRl+>^ 

dx 


( 32 ) 


If we now assume that is known from x and that H, too, is known 
by the unique relation H(H) , equation (52) may be integrated with respect 
to X, and there results, if we introduce in addition the quantity 


and take equation ( 5 I) into consideration 

X^ + f E(H)U^2nR^+n^« 

J X-i 


(33) 


— , .1+n , .5+2n^, ,1+n 
H(x) U(x) R(x) 


(3^) 


©(x) = 
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As a new abbreviation we introduced 


E(H) = 2(1 + n)3(H)H“ (55) 

We plotted this function in figure 6 for the two cases of laminar and 
turbulent flow. 

By way of approximation, we shall now ass\jme constant mean values 
for e(H) and H. This assumption proves correct in a particuleirly 
satisfactory manner for the turbulent case . We then obtain from equa- 
tion ( 5 ^) 


Cl + A r U^^'^R^‘*’“dx ' 

n X, 

^ = ± 

y5t2nj^l*tn 



(36) 


Herein 


A = 


E 


gl+n 


= 2(1 + n)i^ 

H 


(37) 


signifies a mean value suitable for the range < x' < x. 
The integration constant is determined to be 


Cl = \ 

with n = 1/6 if, starting from the point x^, a turbulent boundary 

layer is present. As the laminar or, respectively, turbulent mean value 
for A, we shall choose the value which resvilts when we assume U = 
to be constant and consider the flat plate for fully laminar or fiilly 
turbulent flow. Equation ( 56 ) then becomes, with x^ = 0 and = 0 


0 


P “ 



(59) 


Between the momentxim-loss thickness •flp(l) and the drag coefficient Cf 

of a plate of the length I wetted on one side and approached by a flow 
of the velocity there exists the connection 
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If we put in equation (59) 
tions ( 59 ) and {kO) 


X 


’®p(^) _ 

1 2 

= Z , there follows by comparison of equa- 



A = 



(57a) 


According to the existing flow state, Cf is to be taken for laminar 

or turbulent flow (fig. 7 ) . The notation (57a) offers the additional 
advantage that the surface roughness also can easily be taken into con- 
sideration, merely by substitution of the corresponding Cf values of 

a rough plate. For smooth surfaces, there result the values of table 1 
for the constsint A. 


TABLE 1 .- THE QUAMTITIES n AMD A 
FOR LAMINAR AND TURBULENT FLOW 



Laminar 

(Blasivis) 

Turbulent. 

(Falkner)^ 

n 

1 

1 



6 

A 

O.hhl 

0 .760 X 10 "2 


We now solve equation ( 56 ) for the momentum thickness ■O and 
obtain - if we introduce, in addition, dimensionless quantities anH take 
equation (37^) into consideration - the following expression 



^V. M. Falkner, Aircraft Engineering, I 5 , 19^5 ^ P* 65 . 


Cf 


(hi) 


0.0506 



is valid. 


Ik 
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The integration constant is determined as 


C 




^ 
I I 


1+n 



Summarizing, we repeat once more: If the flow, starting from the initial 

point = 0, is fully laminar or fully turbulent, one has Cx* = 0* 

One has to insert accordingly the laminar or turbulent drag coefficients . 

If a laminar starting length precedes the turbulent boundary layer, the 
value for d in equation (42) is to be taken from the laminar* calcula- 
tion (second formula in equation (42)) . For the laminar case n = 1, 
for the turbulent case n = I /6 is valid. The two-dimensional case 
results if one omits in equations (4l) ajid (42) the radius R. As can 
be seen from equation (4l) , the quantity n does not play a significant role 
in the turbulent case so that one may assume for rough calculations also 
n = 0 . 


(b) The Calculation of the Form Parameter 

We shall determine the form parameter H which is characteristic 
for the velocity profile in the boimdary layer from equation ( 18 ) ; we 
write this equation as follows : 


e(x) ^ = F(ii)r(x) + g(h) 


(45) 


Here © is given according to equation (55) • For the other abbreviations 


and 


r(x) = Q ^ (44) 

U dx 

F(ii) = (H - 1)H 

/ (45) 
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are valid. Fvirthennore, we transform equation (43) ty introducing the 
substitution 


1(B) « = L(H) (1.6) 


and the abbreviation 


K(H) = - = K(L) (47) 

F(H) 


We then obtain a differential equation for the new foim parameter L 

0 (x)H = r(x) - K(L) (48) 


We determined the quantity L by graphical integration of the func- 
tion 1 /f(H) over H. The value L = 0 we have placed in the domain of 
vanishing pressure gradients (constant pressure, flat plate) . In fig- 
ures 8 and 9, Me represent the relation H(L), also for the laminar 
case a(L), and, in figures 10 and 11, the relation K(L) for the cases 
of laminar and turbulent flow. Whereas no influence of the Reynolds nimi- 
ber exists for the function K(L) in the laminar case, that influence 
is rather considerable for the turbulent case . One can show that it is 
possible to give to the equations for calculation of the form parameter 
appearing in the reports of E. Gruschwitz, A. Kehl, and H. C. Garner the 
form of our equation (48)5. The resulting connections for H(L) and K(L) 
have also been plotted in the figures named above . The differences 
between the individual methods are therefore based on the deviations of 
the curves H(L) and K(L) . 

For solution of equation (48), we set up a linear expression for K(L) , 
(compare the figvires 10 and 11) : 


K(L) = a(L - b) (J+9) 

The quantities a and b are obtained, for instance, from table 2. In 
the tvirbulent case, b is, in addition, dependent on the Reynolds nxrni- 
ber and therewith on the length x. 

^A detailed comparison of the individual methods has been carried out 
in an unpublished report of the author. 
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If we substitute equation (^9) irito equation (48), we obtain a 
linear differential equation of the first order for L which we can solve 
in closed form. For this purpose, we meike the substitution 


I 


a 


e 


dx 

0(x) 



u5+2nj^l+n(^\ A 




(50) 


The last relation follows from the fact that in formvila (56) for © the 
denominator is exactly equal to the derivative of the numerator miiltiplied 
by 1/A. The numerical values a/A for smooth svirfaces axe also plotted 
in table 2 . 


Without dealing in detail with the intermediate calculation - for 
which we perform, in addition, a partial integration, in such a manner 
that the derivative dU/d| does no longer appear - we find finally 


L = 


11 Li + in Hlii + 1 
& Ui I 


b(l’) 


In 


U(l*) 


U-, 




(51) 


At the initial point = 0, there is also = 0. The first term in 

equation (5I) then disappears. Especially for the laminar case there 
applies with b = 0 the simple expression 


L = - r in Mih d 11 (52) 

^0 u(|) I 

We shall report later on regarding the initial values for L at 
the stagnation point x^^ = 0 ard at the transition point. 

As can be seen from equation (5l)> 01^6 provide the new vari- 
able I with an arbitrary factor without causing thereby a change in L. 
We ma,y therefore write 


i 




(50a) 
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TABLE 2— THE QUAMTITIES a, b, AMD a/A 
FOR LAMINAR AND TURBULENT FLOW 



Laminar 



L > 0 ^ 

L < 0 ^ 

Turbulent 

a 

2.87 

3.55 

5.04 X 10"^ 

b 

C 

) 

0.071g — - 0.25 

V 

a 

A 

6.5 

8.0 

4.0 


^Pressttre drop. 
^Pressure rise. 


The calculation of | thus will be very simple, since the expression in 
brackets already occxirs in the calculation of the raoment\im-loss thickness 
(eq. (in)). 

A complete calculation is carried out as follows . 


Prescribed; 


General quantities: U(x) , R(x), 

State oj 

^ 1 , Li- 


State of flow n, c^, ^ (laminar or turbulent); initial values Uq^, R^ 


Desired : 

Momentum-loss thickness d(x) ; form parameters L, H(x) . 

First, one calculates the momentum-loss thickness according to equa- 
tion (4l) by performing a simple quadrature. Hence one forms, if one has 
to calculate turbulent flows, the Reynolds number Ud/v with which one 
then calcvilates the qioantity b according to table 2. Next, one deter- 
mines according to equation {^ b .) the new variable | over which one 
integrates the function b - In U/Uq. According to equation (5I), one 

then obtains the form parameter L and, by means of figure 8 or figure 9 
the boundary-layer thickness ratio H. Due to the large exponent a/A 
in equation (50a), the quantity | increases very rapidly with growing 
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length X. This signifies, however, that the term in equa- 

tion (51) loses its significance more and more with increasing distance; 
therewith it is shown that the initial value has, in case of larger 

distances x (for instance in the neighborhood of the separation point), 
only slight influence . The advantage of our method lies in the fact that 
only simple quadratures have to be performed in the individual case. 

Beyond that, no derivatives of the initial values U and R with respect 
to X occur. 


(c) The Initial Values 

The initial values for the momentvun-loss thickness d as well as 
for the form parameter L may be different according to the case to be 
dealt with . 

1. Flow toward a body with stagnation point .- The boxmdary -layer 
calculation is started at the stagnation point = 0 . There U = 0 

and R = 0. From equation (33) then follows immediately 


©0=0 (55) 

In the two-dimensional as well as in the rotationally symmetrical case, 
the potential velocity varies linearly with the length x, that is 


U = cx with c = ( 5 ^) 

The constant c is different for the two-dimensional and for the rota- 
tionally symmetrical case . 

For the variation of the radius, there results also a linearity 
in X, namely 


R = X 


(55) 


The integration constant in equations (3^)^ (36), and (4l) disappears 


Xi = Cl = Cl* = 0 


(56) 
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If we substitute equations ( 5 ^) and (55) into equation (4l), there 
follows, if we put all quantities which additionally enter the formulas 

in the rotationally symmetrical case into braces ^ j^, after a brief 

intermediate calculation 


X — >0 : 


I 


1-n 


1 


>[ 2(2 + n) 



-i(fr 


(57) 


Therein c* = — (— ) signifies the dimensionless expression for c. 

U„\d^/0 

In the laminar case (n = 1) , there results a value different from zero 
for the momentum- loss thickness in the neighborhood of the stagnation 
point 



1 H 

2 \ l 2(3 + {1}) 


(58) 


If one takes into consideration that c* 


write for equation ( 58 ) also 


-i- c and Cf> 

Ux> 


1.328 
— — ^ 



one may 


\f£ % > 1.528 (59) 

In the turbulent case (n = l/6) , the momentum- loss thickness at the 
stagnation point itself (x =0) has the value zero. 


For the momentum- loss thickness ratio at the stagnation point for 
equal velocity increase (dU/dx)Q in the rotationally symmetrical and 

in the two-dimensional case (c = Constant) , one then has 
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1 

. (60) 

c=const J 

If the numerical values for n are substituted, there results in the 
laminar case the value O.867 and in the turbulent case O.816. The exact 
value for the laminar flow which one can calculate from the Hartree pro- 
file is 0.845. 

The initial value of the form parameter L is determined according 
to equation ( 48 ) by putting therein according to equation ( 55 ) ©q = 



Tq = K(Lo) • (61) 

The value Tq is obtained from the fo 3 J.owing boundary-layer determina- 
tion 


0 dU 


© 


To = 11m = lim — = 


►0 ^ X — >0 ^ 2(2 + n) + 1^1 + n^ 

In the laminar case, there resvilts with equation ( 59 ) 


r c «2 


0.220 


3 + 




and from equation ( 49 ) 


It) 


, n.' 


QXL 


3 -. {1} 


(62) 


(62a) 


(65) 


(See table 5.) AlthoiJgh, according to equation (62), in the turbulent 
case a finite value different from zero does result for Fp, it is of no 
^It can be shown that this value follows also from eqviation (52) . 
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help in calculating the value Lq according to equation (6l) since the 
values K(L) for the Reynolds numbers — >0 £ire not known. Since 

V 

the flow at the stagnation point always will be initially laminar, we 
refer to section 5^ c3, where we report on the initial values at the 
transition point. 

For the variation of the form parameter L with the length x at 
the stagnation point, there applies also (as can be derived from equa- 
tion (^8)) 


(S)o = ° 

2 . Flow into a channel . - The boundary layer begins at a point 

= 0 where the velocity, suid likewise the radius in the rotationadly 

symmetrical case, are of finite magnitude. Since again "X-i = = 0 , 

there follows immediately from equations (5^)> ( 58 ), auid (4l) 

= 0 and \ ® (^5) 

and from equation (44), when dU/dx 4 " 

Tk = 0 (66) 

The form parameter then follows from equation (48) as 


K(Lk) = 0 

In the laminar case, this form parameter becomes 


( 6 ?) 


Lk = 0 


( 68 ) 


that is, the value of the flat plate in longitudinal approach flow. 

There is no difference between the two-dimensional and the rotationally 
symmetrical case. (See table 3-) For the turbulent flow, there applies 
what was said above in the discussion of the stagnation-point flow. 

The variation of the form parameter L with x is obtained from 
equation (48) 
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The numerical value a/A is to be taken from table 2. 


( 69 ) 


TABLE 5-- LAMINAR INITIAL VALUES. THE EXACT VALUES - INSOFAR 
THEY DO NOT AGREE WITH THE APPROXIMATION 
VALUES - ARE PUT IN PARENTHESES 



Stagnation-point flow 


Channel flow 


Plane 

Rotat tonally 
symmetrical 


Plane, 

rotat ionally 
symmetrical 

\/F^o 

0.271 
( .292) 

0.235 

( . 247 ) 


0 


.0260 
( .0292) 

.0195 
( .0208) 

Lk 

0 

Ho 

2.25 

(2.22) 

2.52 

(2.5O) 

Hk 

2.60 

°0 

•345 

( .360) 

.520 
( . 324 ) 


.220 


5* Transition laminar - turbulent Behind a certain transition 
region; the laminar boimdary layer is transformed into the txirbxilent layer 
at the point X]_ = x^. From the theory of the origin of turbulence, com- 
pare H. Schlichting (ref. I 7 ), one can find, in agreement with measurements, 
that the transition occurs at places which lie somewhat downstream with 
respect to the velocity maximum. Since the phenomena in the transition 
region have not yet been explored, determination of the initial values 
for ■Q and H or L which are required for the turbulent calculation 
is only approximately possible. What is fundamental in oxir considerations 
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has been shown in figure 12 on the example of the flat plate. Up to the 
point Xy "the boundary layer is laminar and obeys the regularities 

according to Blasius, that is, ^ /jc. The corresponding H-value is 
constant and amounts to H = 2.6o. Starting from the point the 

boundary layer is fully txarbulent . Here applies approximately 

• The corresponding H-value depends, in addition, to 
some extent on the Reynolds nimiber Ito'fl/v prevailing at the transition 
point. According to meastirements , for instance, by F. Schultz-Grunow 
(ref. l8) , also from the similar solutions in case of constant press\ure 
by J. Rotta (ref. 15)^ one finds about 1.2 < H < 1.4, with the H-value 
being the smaller, the larger the Reynolds nimiber. Thus the H-value 
decreases compared to the laminar value . In the transition region, it 
must therefore vary continuously from the laminar to the tiirbulent value. 
In figure 13, we plotted the difference against U^/v. For the saice 

of simplicity, we let Xy^ coincide with x^; then the initial value for 

the momentum- loss thickness of the turbulent calculation is equal to the 
momentum- loss thickness which would result at the point x^ if the flow 

were fully laminar up to this point . We shall therefore put 


^t(xu) = = hM (TO) 

For the H-value, we write 

Ht(xu) = Hu = Hj(xu) - M (71) 

We shall now assume that ^ may be taken from figure 13 also for the 
cases w'ith pressure gradient. This assumption we deem Justified since 
the transition point lies in the proximity of the point of vanishing 
pressure gradient and, thus, the values of the flat plate may be used in 
suitable approximation. Having thus determined Hu, we ascertain according 

to figure 9 the value L^ = L(xu) which then enters equation ( 5 I) as Lq. 

We want here to point out once more that in places lying further downstream 
from transition point the initial value for L is only of slight signif- 
icance; therefore, a somewhat rougher estimate seems Justified when it- is 
a matter of determining the separation point. 


(d) The Separation Point 

Knowledge of the position of the separation point also is important . 
Separation results when the wall-shear stress assmes the value zero. 
Whereas this point in the laminar case corresponding to the prescribed 
velocity profile is fixed by aj = 0 according to equation (l9), it is 
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not yet possible to make a perfectly clear statement regarding the separa- 
ti.on point in the turbulent case . According to the wall- shear- stress 
statements of Ludwieg-Tillmann and Rotta ( compare fig . 5) f "the shear 
stress decreases with increasing value H, however, without attaining 
the value zero. As approximation rule one may, for instance, assume 
according to v. Doenhoff-Tetervin (ref. 11) that the separation starts 
at the earliest when H » 1.8 ”7 and has certainly taken place when H 
has attained the value 2.U. Table 4 presents a compilation of the occxir- 
ring values . 

TABLE 4.- VALUES OF THE FORM PARAMETER 
AT THE SEPARATION POINT 



Laminar 

Turbulent 

H 

4.058 

1.8 to 2.4 

L 

-0 .018 

-0.15 to -0.18 


(e) Examples 

The usefxilness of the method described above is shown on one example 
each of laminar and of turbulent flow. 

1 . Howarth - flow (laminar) .- As example for the laminar flow, we 
choose the well-known Howarth flow (ref. 19)* In this case, the velocity 
distribution is 


iL = 1 - X (72) 

Uoo I 

If equation (72) is substituted into equation (4l) where the radius distri- 
bution R is omitted, the quadrature may be carried out in closed form. 
There results for the momentvim-loss thickness 



^This statement corresponds approximately to Grvischwitz' assumption 
T] = 1 - (u^/u)2 = 0.8. 
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In figure lU, we show the comparison of this approximation with the 
exact values of Howarth. Within drawing accuracy, the agreement is 
perfect . 

For determination of the form parameter, we first calculate the new 
variable | according to equation (50a) 



The form parameter itself also may be calculated analytically according 
to equation (52). Since we find oijrselves, corresponding to the pre- 
scribed velocity distributions, in the region of pressure increase, we 
choose according to table 2 the value a/A =8. We obtain after a brief 
intermediate calculation 


L = i 


0 


i ^ 6 ' d6 * = - -3^ Z'Q 


U d| 


6z 


8J, 


1 - z' 


dZ' 


IT -gg-- 

5j^m+ 8 


(75) 

Taking equation (7^) into consideration, we represented L(x) also in 
figure l4. From these values we determined, with the aid of figiore 8, 
the values a for the wall-shear stress. They, too, are represented in 
figure lU and are there compared with the exact values of Howarth. The 
agreement is quite satisfactory. We also showed the ciorve which results 
according to the Pohlhausen method.® One achieves, as already stated by 
A. Walz, a considerably better determination of the separation point, 
a = 0, if one uses in addition to the momentum theorem also the energy 
theorem, as in our method. 

2. Profile NACA 65(2l6)-222 (turbulent) .- As an example for the 
turbiilent calculation, we choose the profile NACA 65(2l6)-222 (approxi- 
mately) for which measurements (ref. 11) as well as theoretical calcula- 
tions have been carried out. The conditions refer to the upper side of 
the profile placed at an angle of attack a = 10 .1° . The Reynolds num- 
ber is U^Z/v = 2.64 X 10®. The graphical representation of the velocity 
distribution, the momentum-loss thickness, and the form parameter according 
to measurements, likewise the result of our calculation, are given in 


%'his curve we took from A. Walz (ref. 5) • 
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figure 15- We started our calculation at the first point measured, that 
is, x/Z = 0.075- The agreement between our approximation and the meas- 
urement may be called satisfactory .9 For comparison, we also plotted the 
results according to the calculation of v. Doenhoff-Tetervin and according 
to the methods of Garner, Gruschwitz, and Gruschwitz-Kehl . Regarding the 
form parameter H, the last two methods deviate greatly from the other 
methods and from the measvirement . 


9as a. E. V. Doenhoff and N. Tetervin pointed out, the larger differ- 
ences between measurement and calculation for the momentum-loss thickness 
in the neighborhood of the separation point x/Z 0-55 might be caused 
by systematic measuring inaccxrracies . 
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APPENDIX 

DISCUSSION OF THE FORMULA FOR CALCULATION OF 
THE MOMENTUM-LOSS THICKNESS 


As mentioned before, the momentum-loss thickness was calculated, so 
far, from the momentum theorem (16) . With the use of certain simplifica- 
tions, it is possible to find for the calculation of the moment ijm-loss 
thickness from the momentum theorem a formula similar to our equation (56) . 
Two possibilities exist: 

1 . The wall-shear stress is determined according to the laws of the 
flat plate in longitudinal approach flow where the length x is expressed 
by the momentum-loss thickness and the approach-flow velocity Uc 
by the local velocity. For the laminar as well as for the turbulent case, 
one may write 


•’'0 ^ g 

pu2 


If one substitutes, moreover, in the momentum theorem a constant value 
for the boundary- layer thickness ratio H, one can integrate the momentum 
theorem in closed form, as was shown by E. Truckenbrodt in an unpublished 
report (compare H. Schlichting, Boundary-layer theory, page k 30 ; one 
obtains the result 



Cl + (1 + n)a J U®r 1+» dx' 


Ua^l+n 


where a = (l + n)(H + 2 ) - n. 

The corresponding numerical values may be taken from table 5. 
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TABLE 5 




Laminar 

T\arbulent 


n 

1 

i(Blasius) 

k 

i(Falfcner) 

Factor ahead of the 
integral 


0 A 4 l 
A 4 l 
A70 

0 .0160 
.0160 

0 .0076 
.0076 

Exponent of the velocity 
distribution 

3 + 2 n 
a 
b 

5.0 

8.2 

5.0 

4.0 

4.0 

5-53 

5.67 


2 . According to K . Pohlhausen for the laminar case and to E . Buri 
for the tirrbulent case, the wall-shear stress and the boundary- layer ratio 
are functions of the quantity 


0 ^ 
U dx 


that is, 

Jo , £ai£l 

pn^ - 

and 


H = f 2 (r) 


If one puts 


A = 


and 


X* = 



Rl+n 


dx 
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one can find after some intermediate calculations the following equation: 


dA 

dx* 


= 0(r) 


with 


Therein 


A ^ 
U dx* 


0 (r) = (i + n)f3^(r) - ja + n + (i + n)f2(r^r 

represents a universal function for the laminar or, respectively, turbu- 
lent state which, in the first case, may be calculated analytically (for 
instance, for the Pohlhausen polynomials) and, in the second case, may 
be determined from measurements . In both cases , one may find by way of 
approximation a linear connection between $ and F (compare A. Walz, 
Lilienthal-Bericht l4l (l94l) and E. Buri) . 


$ (r) = -bF + c 

If one substitutes this expression in the above equation, the integration 
may be carried out in closed form and yields after a simple intermediate 
calculation 


0 = 



pX 

Cl + c / dsc' 

^^1 

ybj^l+n 


The nxanerical values c and b are also indicated in table 5* (Compare 
H. Schlichting, Boundary- layer theory, pp. 199 and 424.) The constants 
according to equation (56) also have been indicated for comparison. 


Translation by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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Figiire 1.- Survey sketch. 





ZD 3jO 4.0 

Figure 2.- Laminar friction-layer parameters (Hartree profiles). 
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Figure 6. - The function E(H) for laminar and turbulent flow. 
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Figiire 7.- Drag law of the smooth flat plate in longitudinal approach flow. 
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Figures.- Relation between the boundaxy-layer thickness ratio H and 
the wall shear-stress coefficient a, on one hand, and the form 
parameter L on the other for laminar flow. 



Figure 9.- Relation between the boundary-layer thickness ratio H and 
the form parameter L for turbulent flow. 
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Figure 12.- 



Survey sketch from the transition point of the flat plate. 
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I'igure 13.- Variation of the form parameter H in the range of 
transition from laminar to turbulent flow. 
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Figure 15.- Turbulent boundary layer on the profile NACA 65(216)-222 

a = 10.1°, ■— 2.64 X 10®. Measurenaent according to NACA Rep. 772. 
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